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OO , Abstract 

We investigate the structure of constituent quarks and study implications for quark 
On ' distribution functions of hadrons. Constituent quarks are constructed by dressing bare 

p • quarks with Goldstone bosons using the chiral quark model. We calculate resulting 

Qh! corrections to the twist-2 structure functions gi{x) and hi{x). The Goldstone 
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boson fluctuations produce a flavor asymmetry of the quark distribution in the nucleon 
in agreement with experimental data. They also generate significant depolarization ef- 



■ fects which reduce the fraction of the nucleon spin carried by quarks. Corrections to 

the transversity spin structure function hi{x) differ from those to gi{x), and in partic- 
ular we find a large reduction (40%) of the d-quark tensor charge, which is consistent 
with recent lattice calculations. We also study the pion structure function and find 
the momentum fraction carried by the sea quarks in the pion to be considerably larger 
than that in the nucleon. 
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1 Introduction 

High energy experiments such as lepton-hadron deep inelastic scattering and the Drell- 
Yan process provide detailed knowledge of the quark and gluon distributions of hadrons. 
One obtains information about the momentum distribution functions of valence quarks, 
sea quarks and gluons, and their scaling violations are quite consistent with predic- 
tions from perturbative QCD Even under these circumstances, we still have a most 
important problem, namely how to understand the structure functions themselves. 
Perturbative QCD describes their evolution, but cannot predict the distribution 
functions because of their non-perturbative origin. Although lattice QCD studies pro- 
vide some moments of the structure functions p|, a satisfactory understanding is still 
far from being reached. 

On the other hand, much of low-energy hadron phenomenology is quite successfully 
described in terms of the constituent quark picture. At low energies, the sea quark 
and gluon degrees of freedom are assumed to be absorbed into constituent quarks as 
quasi-particles, and hadrons are constructed out of a few such quasi-particles. At the 
high energy scale, sea quarks and gluons reappear and hadrons reveal themselves as 
complex many-body systems with a large number of current quarks and gluons. We 
do not yet have a clear understanding of the connection between constituent quarks at 
low energy and the parton picture at high energy. 

Some attempts have been made to model this connection or at least some aspects 
of it ^, ^ 1^. In such studies, the twist-2 part of a given structure function is 
evaluated within a relativistic quark model, and the resulting structure function is then 
evolved from the low-energy scale, where the constituent quark picture is supposed to 
work, to the high momentum scale with the help of perturbative QCD. Such studies 
may help to approach the deep inelastic scattering data from the point of view of 
low-energy non-perturbative dynamics. For example, the observed large deviations 
of the ratios -F^(x)/-F2 (x) and gi{x)/Fi{x) from the simple parton picture can be 
understood by taking into account the spin-flavor structure of the nucleon [0, Q in 
terms of the structure of constituent quarks. In addition, the difference between the 
M-quark distribution in the pion and the one in the kaon is well reproduced within a 
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chiral model where the constituent quarks play a crucial role. 

Despite its successes, the perturbative QCD evolution becomes problematic when 
it is started from scales far below 1 GeV. In most cases, calculations require that the 
evolution is performed upward from scales as low as /i^ ~ 0.1 — 0.3GeV^ to reproduce 
the high experimental data. In this region the use of perturbative QCD may be 
questionable, even though the difference between results obtained in leading order and 
next-to-leading order from a scale /x^ ~ O.SGeV^ is only about 10%. 

Recent studies of Kulagin et a/. have shown that such difficulties can be overcome 
by taking into account the structure of the constituent quarks. It was found that 
inclusion of the pion dressing and higher mass spectator processes substantially changes 
the shapes of the quark distributions at the given model scale. This modifies the 
normalization of the quark distributions and leads to the correct small- a; behavior by 
introducing Regge exchange. The identification of the results in ref. with the twist- 
2 quark distribution can now be done at a scale of about IGeV^, where the use of 
perturbative QCD for further evolution to high is reasonably justified. 

In the present paper we follow a similar direction, although with different emphasis, 
and extend it to the spin dependent structure functions and to hadrons other than the 
nucleon. At the scale below 1 GeV, the relevant degrees of freedom are assumed to 
be constituent quarks (CQ) and Goldstone (GS) bosons|jlO[. Here, we concentrate 
on the Goldstone boson dressing in the SU{?>)r ® SU{3)l chiral quark model with 
explicit SU (3) flavor symmetry breaking. We study the twist-2 structure functions of 
the nucleon [|l^, namely the unpolarized fi{x), the helicity difference gi{x), and the 
transversity difference hi{x), as well as the pion structure function with inclusion of 
the CQ structure. While Regge high-energy behaviour was found to be important in 
previous studies of the unpolarized structure function, primarily at small x, our main 
focus here will be on the "soft" dynamics mediated by the pseudoscalar meson octet. 
The GS boson cloud influences not only the normalization at the quasi-particle pole of 
the constituent quarks but also changes their spin structure by emitting the GS bosons 
into P-wave states relative to the CQ, thus producing depolarization effects. 



Recently, such a pion dressing has been studied by several authors |T^ in order 
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to account for the violation of the Gottfried sum rule and the nucleon spin structure. 
We first reexamine this previous work systematically. Then we study the chiral-odd 
transversity spin structure function hi{x). Corrections to hi{x) turn out to differ from 
those for gi{x). This difference causes a strong reduction of the d-quark tensor charge, 
while the axial charge remains almost unchanged. This tendency seems to be consistent 
with a recent lattice study [T^ which suggests \6d\ < \Ad\ as opposed to 6u > Au for 
the M-quarks. 

We also show how the pion structure function modifies with inclusion of CQ struc- 
ture, in comparison with the nucleon. We find that the resulting sea quark distribution 
in the pion is substantially enhanced. The contribution to the second moment from 
the sea quarks in the pion is almost twice as large as that in the nucleon. This issue 
will possibly be studied in future experiments. 

This paper is organized as follows. In Section 2, we construct the constituent quark 
Fock states and evaluate their contributions to deep inelastic scattering. We examine 
the Gottfried sum rule and the nucleon spin structure in Section 3, as it was first done 



by Eichten et a/. [12]. In Section 4 we study the chiral-odd transversity spin structure 
function hi{x) and show that the GS boson dressing of quarks changes the simple 
quark model result for hf{x) substantially. Section 5 is focused on the discussion of 
the pion structure function. We emphasize that the inclusion of the GS boson dressing 
naturally leads to an enhancement of the sea quarks in the pion. We draw conclusions 
in the final section with a brief estimate of contributions from multi-pion Fock states. 



2 Structure of constituent quarks 

We start by constructing the constituent quark Fock state using the chiral quark model 
of Manohar and Georgi ||lOl. In this model, constituent quarks couple to the Goldstone 
bosons of spontaneously broken chiral symmetry. The GS bosons, in particular the 
pion, play a crucial role as approximate zero modes of the QCD vacuum. They govern 
the low-energy dynamics at characteristic scales fi^AnfTj ~ IGeV where ~ 93MeV 
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is the pion decay constant. 

Let if) = {u, d, s) be the quark field with Nf = 3 fiavours. The effective interaction 
Lagrangian in leading order is given by 

9A 



Ant = 

with the GS boson matrix field, 
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Here, we set the pseudoscalar decay constant / equal to the pion decay constant /,r and 
start with the quark axial- vector coupling = 1. In their original work, Manohar and 
Georgi took = 0.75 to reproduce the nucleon axial- vector coupling constant within 
the non-relativistic approximation [|l^]. We first adopt = 1 following the large iV^ 



argument from ref. and discuss possible renormalization effects later. 

The d- and s-quarks that enter in eq. are assumed to have already developed 
large dynamical masses as a dynamical consequence of spontaneous chiral symmetry 
breaking. We denote those "bare" but massive states by \uo) and |(io) etc. Once they 
are dressed by GS bosons, we write the constituent u- and (i-quark Fock-states as 
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(3-a) 
(3-b) 



where Z is the renormalization constant for a "bare" constituent quark and |ajp are 
probabilities to find GS bosons in the dressed constituent quark states. The wave 
function renormalization Z only operates on the bare CQ state, because we use the 
renormalized quark- meson coupling constants to calculate a^. In this paper we restrict 
our study to the admixture of one GS boson and truncate the Fock space expansion 
as displayed in eqs. (|3-a] ), ( P-b| ) . We shall give an estimate of two-pion Fock state 
contributions later. 

Within this approximation, the diagrams Figs. 1(a) and 1(b) contribute to the 
structure functions. We use the Infinite Momentum Frame (IMF) to calculate these 





(a) (b) 
Figure 1: 

Diagrams contributing to the constituent quark structure. Fig.l.(a) shows the GS 
boson spectator process, whereas Fig. 1(b) probes the structure of the GS boson itself, 
with the constituent quark being spectator. The virtual photon is depicted by the 
wavy line. Thick sold curves represent the quarks, and the dashed curve the GS boson. 



contributions [jT5|, By virtue of the IMF, the factorization of the subprocess is 

automatic and we may neglect possible off-shell corrections, since all the particles are 
on-mass-shell in this frame. Hence we can use one-dimensional convolution formalism 
throughout the following calculations. Also, working in the IMF removes the so-called 
Z-graph contributions. 

In previous studies we have already obtained the pion dressing corrections to the 
constituent quark distribution function at a low- momentum scale [^, as follows. The 
spin- independent term corresponding to diagram 1(a) is given by 

l3i^)= t^P3c.dy)qJ^] . (4) 

•^^ y \y / 

Here, P{y)ja/i is the splitting function which gives the probability to find a constituent 
quark j carrying the light-cone momentum fraction y together with a spectator GS 
boson {a = vr, K, rj), both of which coming from a parent constituent quark i: 



1 ( 9AmV r 2 {rrij-miyf + k 



2 I U2 

T 



\ f J J y^il-y) [mf-M/J 

where mi,mj,ma are the masses of the i,j- constituent quarks and the pseudoscalar 
meson a, respectively. Mj^ = "^3~^^t _|_ "^^+^t jg ^-j^g invariant mass squared of the final 
state, and m is the average of the constituent quark masses, m = (mj + mj)/2. The 
integral (B|) requires a cutoff which will be specified later. 
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Diagram 1(b) probes the internal structure of the GS bosons. This process gives 
the following contribution: 

qk{x) = f Vk/a (—] Paj/i ( — ) qi{y2) , (6) 



vi 1/2 \yi/ \y2, 

where Paj/i{x) = Pja/i{^ — x). This symmetry relation holds in the IMF if we use 
a momentum cutoff procedure as in ref. [|T^. Here, Vk/a{x) is the quark distribution 
function with flavor k in the GS boson a, with the normalization Jq dxVk/a{x) = 
1. When we calculate the quark distribution explicitly in Section 4 and 5, we will 
use the phenomenological parametrization of the pion structure function at the scale 
IGeV^fl^. As for the kaon and eta, we do not have experimental data and simply use 
the model calculations of ref. |1^ . Ambiguities of the final results arising from the 
choice of the meson structure functions are rather small, at the level of a few %. 
We define the moments of the splitting functions Pja,/i, 

X Pja/i^ — I dx X Pj (^li{x^ . (T) 
^ 

As for the first moments, {Pja/i) = {Paj/i) = {Pa)- In terms of those the renormaliza- 
tion constant Z is then given by 

Z = l-\{P.)-{PK)-\{Pn) . (8) 

We find Z ~ 0.7 using the standard parameter set qa = irt = 350 MeV and a 
cutoff A = 47r/ as specified at the end of this Section. Numerically, about 75% of the 
deviation from Z = 1 comes from the pion dressing and 20 % from the kaons. 

As examples, we explicitly write down the quark distribution functions in the nu- 
cleon using the splitting functions (^). We assume that the bare quark distribution 
functions are given in terms of the valence quark distributions uq{x) and dQ{x). No 
antiquarks are present before the GS boson coupling is turned on. We find 

1 

+ -K/7r ® ® (Mo + do) + Vu/K ® Pku/u ® Uq 

+ \Pu'n/u ®Uq + ^Vu/r, ® P,, ® (uo + do) , (9) 
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d{x) = Zdo{x) + Ptt ® ito + Vd/^ <^P^<^do + -P^ (g) do 

+ ^Vd/n ^Pn^ {Uo + do) + Vd/K ® Pr ® do 

+ ^Pr,®do+^Vd/r,®Pr,®{Uo + do) , (10) 

where the bare distribution functions in the proton have the normahzations 

/ dxuo(x) = 2 , 
Jo 

/ dxdo(x) — 1 . 
Jo 

Here, we use the following short hand notation for the convolution integral: 

P®q^ l"-Piy)q(-] . (11) 

J^ y \yj 

For the antiquaries and the strange quarks we find: 

u{x) = K/tt -Ptt <8) rfo 

+ ®Pn® {Uo + do) + ^K/r, ® Pr, ® {Uq + do) , (12) 



d{x) = Vdi-K ®P-K®uo 



+ ^Vd/n ®Pn® {uo + do) + ^Vd/r, <S> Pr, ® {uo + do) , 



(13) 



s(x) = Pk® (uo + do) + -T^/r, (8) P,, (8) (ito + do) , 

y 



s{x) 



Vs/K ®Pk® {uo + do) + 



{uo + do) 



(14) 
(15) 



We then obtain the valence quark distributions as u^ai = u{x) — u{x), d^aiix) = d{x) — 
d{x), which satisfy the correct normalization with the renormalization constant Z. For 
instance, the first moment of the valence li-quark is 



Jo 



dxUyal {x) 



dx 



Zuo{x) + P^^/d ® Cio + K/tt ® Pjrd/u ® Uo 



+ 2-P«^/« ®uo- K/tt ® P-Ku/d ® do 

1 

+ Vu/K <8) Pku/u ®Uo+ -^Purj/u ® ^0 

2Z + (P,) + 2(P,) + 2^(P,) - (P,) + 2(P;^) + 2^(P,) 



where we have used eq. (||) in the last step. 

Next we evaluate the spin- dependent process and study the spin distribution of the 
nucleon. An analogous calculation as the one leading to eq. (|]) yields: 

AP(y),»/. = (T)' f'"^ K-'".^)'^-'^^ , (16) 

where ^Pja/i = Pjia/i] — Pjia/ii is the difference of probabilities to find helicity +1/2 
quarks minus helicity —1/2 quarks starting from a parent quark with helicity +1/2. 
Note the change of sign in front of /cfi when comparing P and AP. With the standard 
parameter set mentioned above, this spin dependent splitting function gives a negative 
contribution, that is, the helicity-flip process is dominant. In previous studies |T^, 
the relation P = —AP has simply been used, assuming that the GS boson emission 
contributes only to the helicity flip process. However, the ratio of helicity flip and 
non-flip contributions depends on the dynamics. The magnitude of the spin-dependent 
splitting function is rather sensitive to the choice of the cutoff, because the helicity- 
fiip probability is directly proportional to the transverse momentum integral. The GS 
boson dressing given by the formula (|l^) modifies the axial-vector coupling constant qa 
of the nucleon as well as the quark content S of the nucleon spin, as we shall elaborate. 

Now we specify the momentum cutoff function at the quark-GS boson vertex. An 
exponential cutoff is often used in IMF calculations, since such a form factor has the 



correct t and u channel symmetry ||16|| : 



9a 9a exp 



(17) 



4A2 

Here A is the cutoff parameter. This function satisfies the proper symmetry, Pja/i{y) = 
Paj/i{^ — y)- The value of the cutoff is taken to be about A ~ 4:^/^^ ~ IGeV, the 
characteristic scale of spontaneous chiral symmetry breaking. In actual calculations it 
will be determined to reproduce the experimental data of the Gottfried sum rule. 

Using this cutoff function, we show in Fig. 2 the constituent quark-Goldstone boson 
splitting functions, Pja/i{y) and APja/i{y), for the case of a w-quark {i = j = u) 
together with the pion (a = vr) and A = 1.4GeV. For the constituent quark masses, 
we use niu = rud = 360MeV and = 570MeV as typical values guided by NJL model 
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Figure 2: 

Constituent quark-GS boson splitting functions for twist-2 structure functions. 
Pjir/iix), APjT,/i{x) and 6PjT,/i{x) are shown by the solid, dashed and dash-dotted 
curves, respectively. 

calculations and phenomenology. The resulting splitting function Pua/u is peaked 
a.t y < 0.4, and APja/i is negative over most of the y range. If we choose smaller values 
for the cutoff A (say, A = O.SGeV), APja/i becomes positive but stays very small. 
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3 Gottfried sum and nucleon spin in the chiral quark model 
3.1 Update of the Gottfried sum rule 

We start this Section by re-examining the effects of the GS boson fluctuations on 
the Gottfried sum rule and the nucleon spin structure. These were first discussed by 
Eichten et a/.[T2|. The Gottfried sum rule (GSR) is given in terms of the difference of 
the proton and neutron structure functions: 

GSR = \FP{x) - F^{x)] = 1 t dx (u{x) + u{x) - d{x) - d{x)) . (18) 

Jo X 3 Jo ^ ' 

A naive parton model with isospin symmetric sea, u = d, yields GSR = 1/3, whereas 

the experimentally observed value is 0.235±0.026[|l^, which indicates isospin symmetry 

breaking 'u(x) < d{x) in the nucleon sea. If we take into account the shadowing 

correction to extract the neutron structure function from the deuteron data, the GSR 

value is further reduced by about 20-30% With this correction, the empirical 

value of the GSR is reduced to about 0.21 — 0.22. 

We now employ the quark distribution functions of the chiral quark model as given 

in Section 2. Inserting the expressions (|9| JlO|JT^JT3|) into eq. ([18|) , one gets 

GSR = ^i^Z-^-{P^) + {PK) + l{P,)) 

= ^(1-2 (P.)) . (19) 

With the parameters discussed in Section 2, the renormalization constant is found to 
be Z ~ 0.7. Here, we adopt A = 1.4GeV for the cutoff function eq. ([T7|) , and obtain 
Z = 0.67, (P^) = 0.16, {Pk) = 0.09 and (P^) = 0.06. We then find 0.22 for GSR, in 
good agreement with the empirical value. The dominant contribution to the reduction 
of GSR comes from the renormalization of the bare quark state represented by the 
Z factor. The value A = 1.4GeV should be understood as an upper bound of the 
cutoff in the chiral quark model, based on the assumption that the violation of the 
Gottfried sum is entirely given by the Goldstone boson dressing alone. One should 
note, of course, that this A depends on the input constituent quark mass. If we choose 
m ~ 400MeV, the resulting cutoff turns out to be about A ~ 1.1 GeV. 
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3.2 Angular momentum transfer to the meson cloud 

Next we discuss the spin structure of the nucleon. Analysis of all available experimental 
data[^ gives the following decomposition of the nucleon spin in terms of the quark 
spin, with Aq = q ] -q [ +q 1 -q I: 

Am = 0.82 ±0.02, Arf = -0.43 ± 0.02, As = -0.10 ± 0.02, 

S = Au + Arf + As = 0.29 ±0.06 (20) 

at = 5GeV^, to be compared with expectations of the naive quark models (S = 1) or 
with inclusion of relativistic corrections (S ~ 2/3). Using the spin-dependent splitting 
functions (|16|) , the spin fractions of the constituent quarks are modified from their bare 
quark values as follows: 

Au = ZAuo + ^ (AP,) Auo + (AP,) Ado + ^ (AP,) Auo , (21) 
Ad = ZAdo + i (AP,) Ado + (AP.) Auo + \ (AP,) Ado , (22) 

As = {APk) {Auo + Ado) , (23) 

where we have assumed Asq = at the moment. Recall that the first moments 
(AP) of the spin-dependent splitting functions (AP^) are negative, i.e. the spin-fiip 
probability is larger than the spin non-fiip one in the pion emission process. We find 
(AP^) = —0.06, {APk) = —0.01, and (AP^) = —0.03. The pion emission process 
converts part of the spin of the constituent quark into (P-wave) angular momentum of 
the meson cloud. 

In order to obtain numerical results, initial input spin fractions Auo and Ado ^'^^ 
needed. As a first rough estimate, we start from the naive SU{Q) quark model values, 
Auo = |, Ado = which yield the nucleon axial- vector coupling constant Ga = 
Auo — Ado = 5/3 and the nucleon spin S = Auo + Ado = 1 in the absence of the GS 
boson dressing. Inserting these values into eqs. ( pT| , p2| , p3D , we obtain 



Au = 0.86, Ad = -0.29, As = -0.006. 
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With inclusion of the GS boson dressing the nucleon axial-vector coupling constant 
becomes 



Ga = Au- Ad = 1.15 , (24) 

in reasonable agreement with the empirical value. The total quark fraction of the 
nucleon spin is then given by 

^ = Au + Ad + As = 0.56 (25) 



This value is about twice as large as the empirical S = 0.29 ± 0.06 ||2T|. If we allow 
ourselves to vary the momentum space cutoff, it would be possible to obtain a value 
around S ~ 0.3. However, the agreement with the nucleon axial-vector coupling con- 
stant Ga is then lost. We therefore conclude that the nucleon spin problem cannot be 
solved by the GS boson dressing alone, but the depolarization caused by the P-wave 
coupling to the GS bosons is nevertheless a significant effect. 

The value of S, eg. (pSj) , indicates that about 40% of the nucleon spin is carried 
by the relative orbital angular momentum between constituent quarks and their pseu- 
doscalar meson clouds. This result is quite consistent with the recent analysis of the 
nucleon spin decomposition, according to which about ~ 30% of the nucleon spin is 



carried by the quark orbital angular momentum [^. This apparent agreement is of 
some interest. 



3.3 Axial anomaly effects 

Let us now investigate some necessary further steps. In the present framework, we deal 
only with the SU (3) octet of Goldstone bosons to build up the constituent quark struc- 
ture; we do not incorporate the contribution of the axial U{1) anomaly to the nucleon 
spin. After the EMC, SMC and SLAC measurements, many efforts have been made 
to connect the missing nucleon spin with the U{1)a anomaly of QCD. Such anomaly 
contributions to the spin structure of the constituent quarks have been estimated p3| 



using the three flavour Nambu and Jona-Lasinio model which dynamically produces 
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spontaneous chiral symmetry breaking and incorporates the axial U{1) anomaly in 
the form of 't Hooft's effective interaction between quarks Mean field effects from 
such interactions induce a non-trivial spin structure of constituent quarks such that the 
spin fractions change from their "bare" values Auq = 4/3, Ado = — 1/3 and Asq = 
already before GS boson fluctuations are turned on. In a scenario with maximal rj-rj' 
mixing and with inclusion of relativistic bound state wave functions, Yabu et al. find 
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Auo = 0.91, Arfo = -0.37, Asq = -0.12 . (26) 

This result shows the screening of the singlet axial charge induced by the quark- 
antiquark polarization. The dressing with pseudoscalar mesons should be viewed as 
an additional effect beyond the mean field level. It is therefore meaningful not to start 
from the standard SU (6) symmetric relation, but to adopt the values quoted in p6| ) . 
In this case the strange quarks are primordially polarized, and thus the formula ob- 
tained previously is slightly modified. We then find a smaller and reasonable value for 
the nucleon spin as expected. 

S = 0.20, 

Am = 0.60, Ad = -0.28, As = -0.12 . (27) 
However, although the total spin sum S is now consistent with the empirical value 



0.29 ± 0.06 [^1|1) the individual spin fractions of u- and d quarks disagree with the data, 
and hence the resulting nucleon axial-vector coupling constant Ga = 0.85 is about 30% 
smaller than the empirical one. 



3.4 Pion contribution to the cLxial-vector coupling constant 

At this point, we should emphasize that the pion cloud also contributes directly to 
the axial- vector matrix elements illustrated schematically in Fig.3(b), in addition to 
the contributions already discussed above (Fig. 3 (a)) which are nothing but the renor- 
malization and depolarization effects due to GS boson emission. In fact, within chiral 
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models of the nucleon such as the chiral bag model, the contribution from the pion 
cloud can be explicitly calculated[^. In the chiral limit, = 0, the pseudoscalar 



current induced by the pion contributes to the axial-vector matrix element due to the 
zero mass pion pole. One can show p6| that the magnitude of the pion cloud contri- 
bution Fig.3(b) is half of the quark one, Fig.3(a). The nucleon axial- vector coupling 
constant eq. ( p4|) is then modified as 



Gas^ = (P|^7M75r3^|i^) 
= (An - Ad)s^ 

= [{Aua + Aub) - (A4 + Adb)]Sf, 

~ 1.5{Aua - Ada) Sf, , 

where Aqa are the contributions from diagram 3(a), essentially the same as those 
calculated in the previous subsection, and Ag^ comes from the pion pole, Fig.3(b). 
Note that this pion effect contributes selectively only to the isovector axial-vector 

(3) 

matrix elements G\ = Ga- It does not change the flavor singlet axial current and 
therefore leaves the spin S = G^'' untouched, Au^ + Adb = 0. Hence this additional 



effect may well increase Ga to compensate for the small value obtained in eq. (|27|) , with 
the total spin being unchanged. Indeed, we can obtain the following results including 
the pion contribution of Fig. 3(b) in the chiral limit: 

An = 0.93, Ad = -0A3, As = -0.12, 
Ga = 1.35, 
S = 0.21 , 

where we have used eq. ( PUj) as input for the bare spin distribution and introduce the 
cutoff such as to reproduce the GSR. The resulting spin fractions are now reason- 
ably consistent with the empirical ones, eq. (pO]), although the axial- vector coupling is 
somewhat overestimated in the chiral limit, m^r 0. 

In the case of the physical pion with m.^ ~ 140MeV, we do not have the soft pion 
pole, but the axial-vector current itself is changed as 

^m = V'7m75y^ + /.5X- (28) 
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(a) (b) 
Figure 3: 

Contributions to the axial-vector matrix element from the pion cloud. Notations are 
the same as those in Fig.l, and the circle denotes the insertion of the axial- vector 
current. The diagram (a) is basically the same as Fig. 1(a). The diagram (b) is the 
pion cloud effect discussed in the text. 

We must rely on some specific model to calculate the second term of eq. (^), which is 
beyond the purpose of this paper. The results are model dependent , and as a con- 
sequence we cannot draw strong conclusions about the detailed role of such additional 
pion effects. The tendency of the pion cloud to selectively increase the (isovector) axial 
constant Ga is nevertheless obvious. As a sideremark we note that it is also interesting 
to study the Bjorken-x dependence of the quark distribution function arising from such 
pion components, which is expected to be centered in the small-x region. 



3.5 Cutoff dependence 

One might argue that introducing a cutoff as in eq. (17) implies a high degree of 
arbitrariness. On the other hand, a distinction between "soft" and "hard" scales is 
needed, the soft physics being represented by the chiral effective lagrangian. Short 
distance dynamics, not incorporated in the leading order effective interaction (1), must 
then reappear in finite size effects parameterized by form factors with a cutoff related 
to the chiral symmetry breaking scale. For completeness we have investigated other 
forms of the momentum cutoff function in order to estimate uncertainties inherent 
in the cutoff procedure. Consider first a sharp cutoff for the transverse momentum. 
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^(A^ — fcy). To reproduce the observed value of the Gottfried sum the cutoff is the found 
to be A = 0.6GeV. In this case the resulting renormalization factor of the constituent 
quark wave function is Z = 0.63, and the first moments of the splitting functions 
are given by (Pt,) = 0.17 and (APt,) = —0.04. Note that the magnitude of (APt,) is 
smaller than the one obtained with a Gaussian cutoff, that is, the depolarization effect 
by coupling to the pion becomes smaller. When we adopt the input spin fractions of 



Yabu et al, eq. (|2|), we find Au = 0.57, Ad = -0.25, As = -0.11 and S = 0.20, not 
far from the values, eqs. (P?]). 



We have also used a dipole form for the cutoff, [jT^ij^j ■ Here we need A = 
2.4GeV to account for the GSR. The resulting renormalization factor is Z = 0.65. The 
calculated spin fractions carried by quarks are Au = 0.58, Ad = —0.28, As = —0.12 
and S = 0.18. As a common feature of all those different cutoff procedures, we note 
that the leading effect on the Ag's still comes from wave function renormalization 
represented by the Z-factor We obtain similar values for this Z-factor in all cases, 
once the cutoff is fixed to reproduce the GSR. From these studies we learn that the 
cutoff scheme dependence of the results is within 10%. 



3.6 Discussion 

The limitations of the GS boson dressing model to explain the nucleon spin structure 
are obvious. First, no antiquark polarization is produced in this approach, because the 
antiquarks are locked into the GS bosons which carry no spin. However, experimental 
data indicate that the antiquark polarization is not negligible. Several lattice calcula- 
tions also show that the so-called disconnected parts, which may involve the effects of 
the U{1)a anomaly and OZI-violating processes, give substantial and negative contri- 
butions to the quark polarization 0. Such contributions are not incorporated in the 
diagrams of Fig. 1(a) and (b). 

Another restriction is the smallness of the strange quark polarization. In the original 
work by Eichten et al, a large strange quark polarization was obtainedfl^. However, 
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once the SU{3) flavor symmetry breaking is taken into account, the strange quark 
polarization from the GS boson couplings turns out to be small. We recall that the 
numerator of the spin-dependent splitting function is given by (mix— 7712)^— For u, d 
quarks with mi = m2 = rriu, the typical scale of the transverse momentum fc^ is larger 
than the constituent quark mass, and thus the splitting function becomes negative. 
On the other hand, in the strange quark production process, we have mi = and 
7712 = mg. The large strange quark mass cancels the negative contribution from the 
transverse momentum integral. Therefore, the spin-dependent splitting function for 
the strange quark production becomes quite small. 

Of course, GS boson dressing is not the only mechanism in question. We have 
pointed out that it should be seen in combination with other screening effects, involv- 
ing e.g. the axial anomaly and OZI- violating processes, which renormalize the spin 
structure of constituent quarks as in eq. (^6]). The sum of all effects has some analogy 



with the Arima-Horie renormalization of nucleon gf-factors in nuclei|28 



It is nevertheless interesting that the GS bosons carry about 35% of the total nucleon 
spin. Recently, Ji has discussed a gauge invariant decomposition of the nucleon spin 
22| . It is suggested that the nucleon spin is carried in parts by the quark spin (~ 20%), 



the quark (and antiquark) angular momentum (~ 30%), and the gluons (~ 50%) at 
the low energy scale[^. Such a scenario will possibly be tested in future experiments, 
e.g. by deeply virtual Compton scattering off nucleons. 



4 Corrections to the chiral-odd structure function hi{x) 

We have already discussed the spin-independent quark distribution function /i (x) and 
spin-dependent one gi{x) in the previous sections. We now focus on the chiral-odd 



transversity structure function /;,i(a;)|jTT|. The longitudinally polarized structure func- 
tion gi{x), which is chiral-even and can thus be observed in deep inelastic scattering, 
gives the helicity difference of quarks and antiquarks in the longitudinally polarized 
nucleon. On the other hand, the chiral-odd hi{x) structure function corresponds to 
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a target helicity-flip amplitude in the standard helicity basis, and thus does not have 
a simple partonic probability interpretation. The hi{x) structure function provides 
a correlation between the left- and right-handed quarks due to its chiral-odd nature. 
Although hi{x) is not identified with quark spin fraction of the nucleon, the transver- 
sity structure function hi{x) can be understood as the difference of the numbers of 
quarks with eigenvalues -|-1 and —1 of the transverse Pauli-Lubanski operator S't75 
in the transversely polarized nucleon[|lT|. Unlike the axial charge of the gi{x) distri- 
butions, the gluon operators which contribute to hi{x) do not mix quark operators 
under renormalization because of the chirality. The analysis of hi{x) therefore gives 
complementary information on the nucleon spin structure. 

Naive non-relativistic quark models predict fi{x) = gi{x) = hi{x). Relativistic 
quark models such as the MIT bag model give a relation > due to the 

role of the lower components of the quark Dirac wave function ||11||. In such models, the 
following simple expressions are obtained for the first moments of gi{x) and hi{x), the 
axial charge Aq and tensor charge 6q: 



\Aq\ ~ JdMF'-lc'] 
\6q\ ~ Jdh-[F' + ^G'] 



where F and G are the upper and lower components of the quark Dirac spinor. Recent 
lattice simulations indicate that the tensor charge is larger than the axial charge for the 
M-quark, whereas the magnitude of the tensor charge becomes smaller than the axial 
one for the (i-quark results which disagree with simple bag model expectations. 
Here, we study how the pion dressing modifies the hi{x) spin structure function. 

Since the hi{x) structure function is chiral-odd, it cannot be measured by a chiral- 
even probe such as the electro-magnetic interaction. Following the work of loffe and 
KhodjiamirianPUf , we introduce the chiral-odd forward scattering amplitude 



T,{q,p, s)='-j dV"-^ (p, S\T{J,,{i) J(0) + J(0 J,m}\v. S) , (29) 

where J^s = ipjuj^ip is the axial-vector current and J = i/ji/j the scalar current. This 
amplitude is related with the chiral-odd structure function hi{x) through the optical 
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theorem. One can write 

Tf,{q,p, s) = ^s^ - ^^(It^j hi{x)^ (30) 

and obtains the transversity spin distribution function as 

hi[x) = lmhi{x) (31) 

TT 

We use the same techniques as in Section 2 to evaluate the GS boson correction to 
hi{x). The sphtting function for the transversity difference is given by, 

2 

This result tells us that a relation 



P{x) + AP{x) = 26P{x) (33) 
holds among the GS boson corrections, which implies a saturation of Soffer's inequality 



30| in the chiral quark model. The calculated splitting function 5Pja/i is shown in Fig.2 
together with Pja/i and APja/i- Note that 6P{x) is small and positive in the whole 
X region, in contrast to the negative longitudinal splitting function AP{x). The first 
moment of the transversity splitting function is 

{5P^) = 0.05 , 

to be compared with (P^) = 0.16 and (AP^.) = —0.06. The sign difference between 
6P and AP causes a non-negligible modification for the d-quark distribution func- 
tions. Since the numerator of the transversity splitting function is independent of the 
transverse momentum, (SP) is insensitive to the choice of the momentum cutoff. 

In order to estimate the GS boson corrections explicitly, we need the bare quark 
distributions as inputs. Here and for convenience, we use a covariant quark-diquark 
model, which has already been studied in the literaturep], to obtain the input 

distributions for /i(x), gi{x) and hi{x). 

■t-In ref. [D both scalar and axial- vector quark-diquark nucleon vertices are incorporated to calculate 
the unpolarized structure function. Here we take only the scalar vertex for simplicity, although we 
could use the same procedure to calculate the spin-dependent distributions. Hence model parameters 
used here differ from those in ref. Our main concern here is to clarify effects of the Goldstone 
boson dressing on quark distributions. These effects are independent of minor differences in the input 
quark distributions. 
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Assuming 5'[/(4) spin-flavor symmetry, quark distribution functions of bare con- 
stituent quarks can be written as 

Uo{x) = 2Sf{x), doi^x) = Sf{x), 

4 1 
^Uo{x) = -Sg{x), Ado{x) = --Sg{x), 

4 1 
Suo{x) = -Sh{x), 6do{x) = --Sh{x). 

Expressions for Sf{x), Sg{x), Sh{x) are given explicitly in refs. H, The relation 



Sh{x) > Sg{x) is satisfied as in the bag model calculation |jTT|. We do not have to 
elaborate the input distributions here, since the essential features of the following 
results are independent of the shape of these functions. 

As a consequence of the GS boson dressing, the spin- dependent distribution func- 
tions for -u-and (i-quarks are modified in the following way: 

Au{x) = ZAuo{x) + -AP^(g)Auo + AP^®Ado + -AP^(g)Auo 

= ^ZSg + ^AP^®Sg+'^AP^®Sg, (34) 

5uix) = ^ZSh + ^6P^0SH+'^6Pr,0Sh, (35) 
Ad{x) = ZAdo{x) + -AP^ (g) Ado + AP^ (g) Auo + -AP^ (g) Ado , 

= -^ZSg+'^-AP^®Sg-^AP,®Sg, (36) 

Sd{x) = -lz5, + ^5P^®5,--^-(5P,®5, (37) 
6 b 18 

We show in Figs. 4 the spin-dependent u-quark distributions in the proton, Au{x) and 

6u{x), together with the unpolarized distribution u{x). Results with inclusion of GS 

boson fluctuations are given by the solid curves, the ones without the GS bosons by 

the dashed lines. The renormalization of the constituent quark wave function reduces 

these distribution functions from their original distributions, and the small-x region 

is enhanced by the GS bosons. However, the relative magnitudes of u{x), Au{x) and 

6u{x) are not modified very much. The first moments of Au{x) and 6u{x) are shown 

in Table 1 with corresponding moments of the original distribution functions. 
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hi{x) and gi{x) are similar at a low starting scale ~ O.SGeV as suggested by simple 
quark models, hi{x) becomes much smaller than gi{x) for small x < 0.1 at the experi- 
mentally accessible scale due to the difference in the evolution. On the other hand, 
the GS boson cloud effects studied in the present paper show up at intermediate values 
of the Bjorken variable, 0.1 < a; < 0.5. In particular the rf-quark contribution to hi{x) 
is much reduced around x = 0.2 ~ 0.3 from the simple quark model estimate. The GS 
boson effect also reduces the ci-quark tensor charge very much as already pointed out. 
We conclude that the GS boson cloud effects, if existent, should be seen experimentally 
in the intermediate x region. 

In this paper we omit the Regge-type contributions at small-x, although significant 
effects for the unpolarized structure function were found in ref. . It is clearly an inter- 
esting point, still under investigation, to combine the GS boson cloud effects presented 
here and Regge phenomenology in order to reach a unified picture of non-perturbative 
features at all x. Small-x behaviour was studied for gi (x) and hi (x) [0 within the 
leading logarithmic approximation of perturbative QCD, and interesting differences 
were found. 



5 Sea quark distributions in the pion and the nucleon 

In the framework of the chiral quark model, constituent quarks and Goldstone bosons 
are treated as quasi-particles. They appear structure-less at the tree-level. As shown 
in previous sections, however, the dressing of quarks plays an important role in under- 
standing the nucleon structure. Here we try to describe the pion in the similar way, 
and see how the pion structure function changes when including higher order correc- 
tions derived from the effective Lagrangian of the chiral quark model. We specially 
concentrate on the sea quarks in the nucleon and the pion and show that the sea quark 
momentum fraction in the pion is expected to be significantly larger than that of the 
nucleon. 

Typical diagrams contributing to the pion-virtual photon forward scattering am- 
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Figure 6: 

Typical diagrams contributing to the substructure of the pion in the chiral quark model 

plitude are displayed in Fig. 6. The first diagram, Fig.6(a), represents the elementary 
pion process {i.e. elastic Compton scattering on the pion) and does not contribute to 
the leading twist distribution. The second graph, Fig. 6(b), can be calculated within 
the CQM as a first non-trivial contribution. We denote this distribution function by 
Vq/Tr{x) and associate it with the valence quark and antiquark content of the pion. The 
third and forth diagrams, Fig.6(c)-(d), correspond to the case in which the quark or 
the antiquark in the second diagram is dressed by the Goldstone bosons. These con- 
tributions can be evaluated by using the procedure developed in the Section 2. Such 
processes determine the sea quark distribution of the pion in the chiral quark model. 

We have already given expressions for the quark distribution of the nucleon in eqs. 
(^,|TD|,[T^,|T3p . The sea quark distribution in the nucleon is 



Sn{x) = 2U{x) + 2d{x) + s{x) + s{x) . (38) 

Such sea quark distribution functions arise from the diagram. Fig. 1(b), in which the 
constituent quarks act as spectators. The contribution from this diagram is soft and 
concentrated in the small x region as a result of the double convolution. 

On the other hand, for the vr'^ case in which uq and are bare quark distributions. 
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the non-valence components that arise from the dressing are given by 

u{x) = P^^do + ® P. ® K + do) + ® P„ ® (% + do) , (39) 

d{x) = Pn uo + ^Vd/n ® {uo + do) + ^Vd/^ ® {uo + do) , (40) 

4 

s{x) = Pk®uo + Vs/k ®PK®do + -Vs/n ®Pr,® {uo + do) , 

y 

4 

S{x) = PK®do + Vs/K ®Pk®Uo+ -K/r, ® Pr, ® {uo + ^o) • (41) 

y 

For the bare quark distribution we have Uo{x) = do{x) = Vg/^(x) by definition. Thus, 
one can write the sea quark distribution in the positively charged pion as 

S^{x) = 2u{x) + 2d{x) + s{x) + s{x) . (42) 

Note that both diagrams 1(a) and 1(b) contribute to the sea quark distributions in the 
meson case because of its qq structure. The diagram, Fig. 1(a), gives harder distribution 
functions. This leads to an interesting difference between the nucleon and the pion. 

Let us consider the momentum fraction carried by the sea quarks. From eqs.(|T^, 
^) and (^,1^,0), the second moments of the sea quark distributions in the nucleon 



and the pion are expressed as 

(xSn) = 2{xVg/^) ^{xP^ j /i) {x{uo + do))N H , (43) 

{xS^) = 2{xP,^/i){2xVg/^)^ + ■ ■ ■ , (44) 

where we only write down the dominant contributions. Here (2xVq/7r)7r and {x{uo + 
do))N are the second moments of the bare valence distributions in the pion and nucleon, 
respectively. 

One can establish a simple relation between the GS-boson-induced sea quark dis- 
tributions in the nucleon and the pion. Neglecting the non-leading terms in (^) and 
which are about 20% of the first terms, and taking a ratio, we obtain 

{xS^) 2{xPj^/i) 



{xSn) {xPnj/i) {x{Uo + do))N 

^ , (45) 

{x{uo + do))N 
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where we use the CQM parameters determined in Section 3. This approximate relation 
indicates that the sea quark momentum fraction in the pion is substantially larger than 
that in the nucleon, because {x{uo + c?o))7v must be less than 1. 

Here we simply take {x{uo + do))^ ~ 0.6 from the parametrization of the nucleon 
structure function ||3^ at the scale = IGeV^, at which this relation is assumed to 



be meaningful. Note that 5*,^ and Sn are the sea quark distributions generated by the 
GS boson dressing alone, and cannot be simply identified with the sea quarks of the 
standard parametrization which may also include sea partons from the perturbative 
gluon radiation. In ref. the sea quarks in the nucleon carry the momentum fraction 



~ 0.15 at IGeV^. As an example, suppose that half of this value is generated by the 
GS boson dressing, i.e. {xSn) ~ 0.7. Then we get from eq. (0), 

{xS^) ~ 1.7 , 

a large contribution. The total sea quark momentum fraction in the pion combined 
with the sea quarks generated by the gluon radiation may be more than 0.2. Such a 
large sea quark distribution cannot be generated perturbatively at ~ IGeV^. The 
systematic analysis of the sea quarks in the pion might be a key to understand the 
non-perturbative structure of the constituent quarks. 

Even if sea partons exist in the nucleon and/or pion before the GS boson dressing, 
the relation (45) remains qualitatively unchanged. For instance, in such a case, eq. (^) 
is rewritten as 

{xSn) = 2{xVg/j,)^{xP^j/i){x{uo + do))N + 3{xPj^/i){x{uo + do))N H • 

The contribution from the second term is estimated to be 20% at most of the first 



term according to the parametrization of the antiquark distribution P9||, although the 
second term arises from the process of Fig. 1(a). 

At present, the sea quark distribution in the pion is not determined experimentally. 
In the analysis of ref. [Q, both small (10%) and large (20%) sea quark distributions in 
the pion lead to an equally good fit to the Drell-Yan data. Future experiments should 
clarify this uncertainty. A precise determination of the pion sea quark distribution may 
be possible by using the pion-deuteron Drell-Yan experiment] 
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6 Summary and discussions 

We have studied the Goldstone boson dressing of constituent quarks in the chiral SU (3) 
quark model and its implications for nucleon and pion structure functions. The GS bo- 
son contributions, in which the pion is dominant, produce substantial renormalization 
effects of the constituent quark properties. Such a renormalization yields the asym- 
metry of the u and d quark distribution in the nucleon. The resulting value for the 
Gottfried sum is consistent with the observed data. Similar results are obtained in 
hadronic model calculations with the meson clouds []TB[. 

The CQ spin structure is also directly modified by the GS boson fluctuations. 
They change the CQ spin structure by emitting the GS boson in a P-wave relative 
to the quark. But here the situation is more subtle. The spin-dependent splitting 
function AP{x) is sensitive to the value of the momentum cutoff used in the chiral 
quark model. We have fixed this cutoff to reproduce the empirical Gottfried sum, and 
found that the spin fractions of the dressed CQ are reduced by about 10-30% from the 
ones of the undressed CQ. This depolarization effect is a consequence of spin being 
converted to P-wave orbital angular momentum of the meson cloud surrounding the 
quarks. It is difficult, however, to decrease the spin fraction further, when we keep 
simultaneously the nucleon axial-vector coupling constant Ga at its measured value. 
Hence we conclude that the nucleon spin problem is not solved by the GS boson dressing 
alone. On the other hand, the calculated depolarization in the present study (30-40%) is 
quite consistent with the recent decomposition of the nucleon spin^^ when identified 



with the contribution from the quark orbital angular momentum. When combined 
with other contributions, namely the axial U{1) anomaly and the meson cloud effects 
illustrated in Fig. 3(b), to be identified with the disconnected contributions on the 
lattice QCD simulations, the present results go altogether into the right direction of 
providing a reasonable description of the nucleon spin structure. 

We briefly comment on additional corrections from the multi-pion cloud of the 
constituent quark. Consider the following two-pion Fock space component of the con- 
stituent M- quark: 
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Figure 7: 

Typical diagrams contributing to the constituent quark structure function with the two 
GS bosons. Notations are the same as those of Fig.l. 

Such contributions are described by the diagram, Fig. 7. Here we incorporate the 
ladder diagrams and neglect the crossed pion loop diagrams, which give much softer 
contributions to the structure function T5[ . 

Due to factorization in the IMF, corrections from such diagrams are expressed by 
the already known splitting functions in most cases. We find the following two-pion 
contribution to the renormalization constant: 



= 0.03 



(47) 



which is about 10% of the total single GS boson contribution. We conclude that the 
inclusion of multi-pion Fock states is not expected to change our results significantly. 
The perturbative treatment of the GS boson cloud should therefore be justified. 

We have also evaluated the GS boson correction to the chiral-odd hi{x) structure 
function. In comparison with the result for gi{x), where the CQ-GS boson splitting 
function AP is negative, we find a positive contribution to the hi{x) structure function. 
With inclusion of the renormalization effect, the transversity spin structure function 
for the (i-quark is shown to be reduced significantly. We have found the relation 
\6d\ < \ Ad\, which seems to agree with recent lattice QCD results p!3 . 



One may wonder how the QCD evolution of the tensor charge influences this dis- 



cussion. In fact, the tensor charge decreases by the QCD evolution [^, However, 
even if one evolves from the very low scale fi"^ = 0.2GeV^ to a few GeV^, the reduction 
of the tensor charge is at most 10%. We may therefore expect that the evolution does 
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not change the relative magnitude of Ag and 5q very much. We need some other dy- 
namical mechanism to get An < 5u and |A(i| > \5d\, such as the GS boson dressing of 
the constituent quark, if the lattice result is confirmed. 

Another interesting result has been obtained for the pion structure function. Sea 
quarks in the pion appear by direct emission of the GS bosons, while the nucleon 
case requires the second order process illustrated in Fig. 1(b). This picture naturally 
gives a substantial enhancement of the sea quark distribution function in the pion. 
Such an enhancement can be checked in a future Drell-Yan experiment |[40|| . Our model 
calculations are done at the low energy scale ~ IGeV^, where the chiral quark model is 
supposed to make sense. To make a prediction for forthcoming experiments, we must 
carry out the evolution from the low energy scale. However, the singlet evolution 
equation needs the gluon momentum fraction in the pion at the model scale, which we 
cannot estimate reliably at the moment. 

We have studied aspects of constituent quark structure in deep inelastic processes. 
Non-perturbative features of the constituent quark, i.e. Goldstone boson fluctuations 
around the CQ, lead to screening effects and to characteristic changes of the CQ spin 
structure. Experimental tests of the nucleon transversity spin distribution function and 
the pion sea quark distribution will provide new insights into the role of the constituent 
quarks as quasi-particles. 

In closing we mention that, recently, considerable interest has been focused on the 
possible experimental observation of such a "soft" pion cloud surrounding the nucleon. 
Semi-inclusive deep inelastic scattering processes with pions, nucleons or A-isobar in 
the final state are studied as promising options by several authors [p]] . 
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